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Abst ract 

A case of the matrix Kummer relation of Herz ( 19551 ) based on the Pearson VII type 
matrix model is derived in this paper. As a consequence, the polynomial Pearson 
VII configuration density is obtained and this set the corresponding exact inference 
as a solvable aspect in shape theory. 



Introduction 



Start assuming that a given sample of n "figures" , comprised in N landmarks (" anatomical" 
points) in K dimensions, and summarized in an N x K matrix X, belongs to certain matrix 
variate distribution with unknown scale and location parameters. The statistical theory of 
shape pursues the distribution of the transforming X after filtering out some non important 
geometrical aspect of the original figures, such as the scale, the position, the rotation, 
the uniform shear, and so on; then the so called shape of the population object should 
be inferred and some comparisons among population shapes could be performed, among 
many others statistical comparisons in the shape space, instead of comparisons in the very 
noised non transformed Euclidean space. As usual, the classical theory of shape, with 
different ge ometrical filters, as the Eucli dean or affine, was based on Gaussian samples, see 
for example iGoodall and Mardia I ([19931 ) and the references therein. Then, the non normal 
applications demanded general assumptions for the samples, and the so termed generalized 
shape t heory was set under elliptical models. U nder the Euclidean filter we can menti on the 
work of iDiaz-Garci'a and Caro-Lopera ( 2010l ). Di'az-Garci'a and Caro-Lopera ( 2012f ): they 
find the density of all geometrical information about the elliptical random X which remains 
after removing the scale, the position and the rotation. Finally, if an affine filter is applied, 
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in order to remove from X all geometrical information of scale, position and uniform shear, 
Caro-Lopera et al (|2010h obtained the so term configuration density of X. All the above 
densities are expanded in terms of the well known zona l polynomials, in a series of papers 
by A.T. James in 60's, see for example iMuirheadl (|l982h . 

The transition of the Gaussian shape theory to the elliptical shape theory demanded 
some advances in integration involving zonal polynomials (see for example lCaro-Lopera et al 
jut important problems remain, the computability of the shape de nsities 



In t his paper we focus in alternative s for su ch p roblems. It is easy to see i n Goodall an d Mardia 



(1993), Di az-Garci'a and Caro-Lopera ( 2010() and lDfaz-Garcfa and Caro-Lopera ( 2012 ) that 
the structure of shape densities under Euclidean transformations involves series of zonal 
polynomials whi ch heritages the difficulties f or computations of the classical hypergeometric 



series studied by IKoev and Edelman ( 20061). However, a clas s of th e generalized confluent 

3S ofK - 



type series of the configuration densities of ICaro-Lopera et al ( 2010t ) can be handled in or- 
der to transformed the series into polynomials, and then the addressed open problems for 
computations of the series can be avoided. 

The configuration density under Kotz type samples (including Gaussian) has this prop- 
erty, and the infe r ence c an be performed with polynomials instead of infinite series, see 
Caro-Lopera et al ( 20091). The source for this property resides in a generalization of the 



Kummcr relation of iHerz I (|1955l ) 



This motivates the present work, claiming that there is a similar Kummcr relation based 
on a Pearson VII distribution, which under certain restriction of the parameters in the 
associated Pearson VII configuration density, it can be turned in a polynomial density. Then 
the inference can be performed easily by working with the exact likelihood which is written in 
terms of very low zonal polynomials. In fact the exact densities can be write down by using 
formulae for those polynomials; for example, in the planar shape theory (the most classical 
ap plications resides in the s tudy of figures in Sft 2 ), we can use directl y the formulae given 
by Caro-Lopera et al ( 20071) instead of the numerical approaches by IKoev and Edelman 



2006) in order to perform some analytical properties of the exact density. 



This discussion is placed in this paper as follows: section [2] defines a Pearson VII type 
series and finds an integral representation which lead to a matrix Kummer type relation 
which we call Kummer -Pearson VII relation; then by applying some general properties 
studied by IHerz (1955), the equality is extended for the required domains in the shape 
theory context. Finally, Kummer-Pearson VII relation gives the finiteness of the Pearson 
VII configuration density in section [3] 



2 Matrix Kummer-Pearson VII relation 



Recall that the matrix Kummer relation (due to IHerz I (|1955l )) states that 
iFi(a;c;X) = ctr(X) \Fi{c — a; c; — X) and; 



(1) 



see also IMuirheadl (|1982[) . 

Now, let X > be an m x to positive definite matrix, then define 

/(t,X)^(a), 



L P x (/(t, X) : a; c; X) = £ l^f> £ f C T (X), 



t=o 



(2) 



where the function fit, X) is independent of r, r = {t\, ■ ■ ■ ,t m ), t\ > t% ■ ■ ■ > t m > 0, is a 
partition of t, 



(/9)r = n(/8-5(*-l) 
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and 

(b) t = 6(6 + l)---(6 + i-l), (6)o = l. 

Then, using this notation we see that the Kummer relation ([T]) is a particular case of a 
general type of expressions with the following form 

1P1 (f(t, X) : a; c; X) = wpC)^ (<?(«, X) : c - a; c; h(X)) , (3) 

where the functions v,g and h are uniquely determined by the particular function / and 
according to the domain of the parameters a and c and the matrix X. 

First, we consider an integral representation of the left hand side of ([3]) under the model 
/(f,X) = (6) t . 

Theorem 2.1. Let X < I, Re(a) > \{m— 1), Re(c) > |(m — 1) and Re(c-a) > \{m — 1). 
T/ien /or suitable reals b and d, we have that 

b-t . „. „. v\ _ r m (c) 



l P 1 ((b) t d- b - t :a;c;X 



r m (a)r m (c - a) 

/ (d-tr(XY)) _f, |Y| a -( m+1) / 2 |I- Y| c - Q - (m+1 )/ 2 (dY). (4) 



/0<Y<I 

Proof. First, we use a zonal polynomial expansion 

(d-tr(XY))" 6 = ^W-i_[tr(XY)] 



t! 

t=o 



-b-t 



t\ 



Then integrating term by term using iMuirheadl (| 19821 theorem 7.2.10), we have that 
(d - tr(XY))~ b |Y| a -( m+1 )/ 2 |I - Y| c - Q -(" l+1 )/ 2 (oY) 

V ^ td ~'" t V / |Y| a -( m+1 )/ 2 |I-Y| c - a -( ro+1 )/ 2 C r (XY)(dY) 



/ 

00 ,,. ,_ b _ t 



_ {b) t d t ^ (a) r r m (a)r m (c - a) 

- h r ^ a(x) 

_ r m (a)r m (c-a) x - {b) t d- b - t ^ (a) T 

= rm(a ^ T " ( r a) 1^1 ((&)**-*-* = «; x) , 
r m (c) 

and the required result follows. □ 

Now, we derive the version of ([1]) but based on a Pearson VII type model, we call this 
expression, Kummer-Pearson VII relation. 

Theorem 2.2. Let X > 0, Re(a) > ±(m-l), Re(c) > |(m- 1) and iJe(c-a) > |(m-l). 
T/ien /or suitable reals b and d, the Kummer-Pearson VII relation is given by 

1 F 1 ((6) t d- fc -* :a;c;X) = (d-trX)- b 1 P 1 ((6) t (d - tr X)"* : c - a; c; -X) . (5) 
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Proof. Consider W = I — Y in then we obtain 
1 P 1 ((b) t d- b ~ t :a;c;X 



/. / v , _ r m (c) 



X 



T m (a)T m (c - a) 

f (d- tr[X(I - W)]) -6 |w| c - a -< m+1 )/ 2 |I - W| a - (m+1 ^ 2 ((iW) 



/0<W<I 

r m (c) 



X 



r m (a)r m (c - a) 

/ (d- trX - tr(-XW))~ b | W | c - a -( m+1 )/ 2 |I - W\ a - {m+1 ^ 2 (dW) 

•/0<W<I™ 



'0<W<I 

_ r m (c) 
r m (a)r m (c - a) 

r m (a)T m (c-a) f v\-b-t 

l "l (o)t (a - tr X) : c - a; c; 



r m (c) 

which is the required result. □ 

The rea der can comp are theorem l 2 . 2 1 ( and its proof) wi th the Kummer relation (and its proof 
given by IHerz So the analysis of IHerz I ()l955l ) for extending the above relations for 



other values of the parameters, holds in the Kummer-Pearson VII relation too. 

Explicitly, we proved that the integral representation Q of iPi ((b) t ■ a; c; X) holds for 
Re(X) < I (by analytic continuation), Re(a) > \(m — 1), Re(c) > \(rn — 1) and Re(c-a) > 
\(m — 1). Then by a suitable modification of the arguments in IHerz I ( 19551 ). we can extend 
the domain of as follows 

Theorem 2.3. Re(X) > 0, Re(a) > ~(m — 1) and Re(c) > |(m — 1). Then for suitable 
complex numbers b and d, the Kummer-Pearson VII relation is given by 

xPi ((b) t d- b - t : a; c; X) = (d - tr X)" b 1 P 1 ((&)« (d - tr X)~* : c - a; c; -x) . (6) 

The above relation is important in shape theory applications, in the so called polynomial 
Pearson VII configuration density. 

3 Polynomial Pearson VII Configuration Density 

Our motivation for studying finite shape densities, comes from the computations of hyperge- 
ometric series type involved i n these distributions . It is k nown, that the zonal polynomials 
are computable very fast by iKoev and Edelman ( 20061 ). but the problem now resides in 



the converg ence and the truncat i on of the series of zonal polynomials. In fact, in the same 
reference of Koev and Edelman (200G) we read: 



"Several problems remain open, among them automatic detection of convergence .... and 
it is unclear how to tell when convergence sets in. Another open problem is to determine 
the best way to truncate the series. " 

Thus the implicit numerical difficulties for truncation of any configuration d ensity moti- 



vate tw o areas of investigation: first, continue the numerical approach started by (jKoev and Edelman 



(2006)) with the confluent hypergeometric functions and extend it to the case of some con- 
figuration series type Kotz, Pearson VII, Bessel, Logistic, for example; or second , propo se a 
theoretical approach for solving the problem analytically fsee lCaro-Lopera et al (2009])). 



We study now the second question corresponding to the polynomial Pearson VII config- 
uration density. 
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Recall that a p x n random matrix X is said to have a matrix variate symmetric Pearson 
type VII distribution with parameters s,R € $1, M : p X n, £ : p x p, $ : n x n with R > 0, 
s > np/2, £ > 0, and $ > if its probability density function is 

r(s) A , trtx-MyE-Hx-M)*- 1 



(7ri?)"P/ 2 r (s - 2j£) |53|"/2|$|p/2 V i? 

When s — (np + R)/2, X is said to have a matrix variate ^-distribution with R degrees 
of freedom. A nd in this case, if R = 1 , then X is said to have a matrix variate Cauchy 
distribution, sejCaro-Lopera et aA (l2010h ■ 

Then, by Caro-Lopera et aA (|2010h we have that (see Goodall and Mardia ( 19931 ) for 
the gaussian case), 

Lemma 3.1. Let be 

= MS! = N ~ 1 

X=±U'V- l w"Z- 1 U(U"£- l U)-\ c=£, (8) 

K(N-l) tr(^'S-V) 
b = s-^ '-, d=l+ ^ R (9) 

If Y ~ En-ixk(pn-ixk i ^n-ixN-i®Iki h), for S > 0, then the non-isotropic noncentral 
Pearson type VII configuration density is given by 

A 1 P 1 ((b) t d- b - t :a;c;X), (10) 

where i-Pi(-) has been defined in 0). 

Unfortunately, the above configuration density with general form A i -Pi (/(f) : a; c; X) is 
an infinite series, given that a = and c = -f- are positive (recall that N is the number of 
landmarks, K is de dimension and N — K — 1 > 1). So a truncat i on is needed for performing 
inference when the modified algorithms of iKoev and Edelman I (|2006l ) are used. 

But the above series can be turned into a polynomial if we use the following basic principle 
of lCaro-Lopera et 3 (l2009h . 

Lemma 3.2. Let N — K — I > 1 as usual, and consider the definition of l-Pi(-) in (0). The 
infinite configuration density has the general form 

CD 1= w(K,N,X) iPj (/(f,X) : ^-1;^;X 

for suitable functions: w(-), independent oft and r but dependent of K,N and X; and /(•)> 
independent of r, but dependent of t and possibly of X (it depends on the generator elliptical 
function, compare with the particular Pearson VII case of lemma \3.1\) . Then, according 
to (0), if the dimension K is even (odd) and the number of landmarks N is odd (even), 
respectively, then the equivalent configuration density 

( ( N - 1 K\ K 

CD 2 = w{K, N, XMX) 1 P l [g(t, X) : - (^^— - - J ; - ; h(X) 

is a polynomial of degree K ( jV 2 ~ 1 — -y) m f/ie latent roots of the matrix X (otherwise the 
series is infinite); where v,g and h are functions understood in the context of and depends 
on the elliptical generator function. 
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Given an elliptical configuration density CD\ indexed by the function /(•) and based 
in the fact that a = > 0,c = If > 0, the crucial point here consists of finding an 

integral representation valid for c — a = — N ~^~ l < 0, which will lead to an equivalent 
elliptical configuration density CD2 indexed by some function <?(•). Then the finiteness of 
CD-2 follows from K even (odd) and N odd (even), respectively. 

In particular, for the Pearson VII generator function, the referred polynomial density is 
provided by applying the new Kummer-Pearson VII relation of theorem l2.3[ via lemma [3~2l 
in lemma [3TTI 

Theorem 3.1. Let be A, a, b, c, d and X defined by (7p-(G|). 

If Y ~ En~ixk (hn-ixk i 'En-ixN-i <8> I_rt, h), X > 0, K is even ( odd) and N is odd 
(even), respectively, then the polynomial non-isotropic noncentral Pearson type VII configu- 
ration density is given by 

(d-tTXy b 1 P 1 f(&) t (d-trX)~* : c-a;c;-X) (11) 

and it is a polynomial of degree K ( N Z l — W) in the latent roots o/X. 

Proof. The proof is trivial, just start with the infinite configuration density (|10|| : 

A 1 P 1 ((b) t d- b - t :a;c;X), 

where A, a,b, c, d and X are given by ©-(JS]). Then apply ^ and the result follows. Note 
that the finiteness follows from Lemma (|3.2[) by noting that c — a — — Af ~^ r ~ 1 is a negative 
integer, when K is even (odd) and N is odd (even). □ 



The principle of lemma 13.21 is based on a known property of the hypergeometric series 
easily extended to series of the type \P\{f(t, X) : a, c; X), see ([2]), which states that if a is a 
negative integer or a negative half integer, the series vanishes in a polynomial. Then if wc 
have an application following a confluent distribution type, it is a polynomial, always that 
the parameter a accepts the addressed special domain, otherwise the distribution is a series 
of zonal polynomials and the open problems for its computability appears. This last case 
occurs for example in the general configuration density, which is an infinite confluent series 
type because a — Z 1 > 0, then it is not a trivial to turn that series into polynomials, 
because the application docs not allow a negative parameter. Then the lemma l3~2l gives a 
special subclass by selecting the number of landmarks (even, odd) given the dimension (odd, 
even), otherwise the configuration density remains a series, then the associated Rummer 
relation must be obtained in order to transform the numerator parameter a into c — a 
which is a negative half integer as required. So, in the context of shape theory under affine 
transformations, the study of Kummer relation type p lays an impo rtant role. It is easy to 
check that the classical Kummer relati o n first derive by Irlerz (1955), and set in the context 
of zonal polynomials by IConstantine I (1963) is related with a Gaussian kernel, but if wc 
want to obtain some non Gaussian polynomial densities (under the explained restrictions of 
N and K), then we need to derive new Kummer relations, it was the case of certain class of 
Kotz configuration densities (which includes the classical Gaus s ian), it need the derivation 
of the associated Kummer Kotz relation, see Caro-Lopera et al ( 20091 ). Thus, in the case of 



the Pearson polynomial configuration density the corresponding Kummer relation was the 
key point for transforming series. 

The above discussion opens related problems in some special topics of matrix variate 
analysis involving confluent matrix with special domains for the parameters, or general 
hypergeometric series type with more or equal Pochhammer symbols in the numerator, 
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several examples of this situations, which demands new developments for Euler re l ations 
and similar ones, can be inferred from some distributions proposed in iMuirheadl (|1982 , 
chapters 8-11). 
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